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[i] Deuterium excess (d) is interpreted in conventional paleoclimate reconstructions as a 
tracer of oceanic source region conditions, such as temperature, where precipitation 
originates. Previous studies have adopted coisotopic approaches (using both 5 ls O and d) to 
estimate past changes in both site and oceanic source temperatures for ice core sites using 
empirical relationships derived from conceptual distillation models, particularly Mixed 
Cloud Isotopic Models (MCIMs). However, the relationship between d and oceanic surface 
conditions remains unclear in past contexts. We investigate this climate-isotope 
relationship for sites in Greenland and Antarctica using multiple simulations of the water 
isotope-enabled Goddard Institute for Space Studies ModelE-R general circulation model 
and apply a novel suite of model vapor source distribution (VSD) tracers to assess d as a 
proxy for source temperature variability under a range of climatic conditions. Simulated 
average source temperatures determined by the VSDs are compared to synthetic source 
temperature estimates calculated using MCIM equations linking d to source region 
conditions. We show that although deuterium excess is generally a faithful tracer of source 
temperatures as estimated by the MCIM approach, large discrepancies in the isotope-climate 
relationship occur around Greenland during the Last Glacial Maximum simulation, when 
precipitation seasonality and moisture source regions were notably different from the present. 

This identified sensitivity in d as a source temperature proxy suggests that quantitative climate 
reconstructions from deuterium excess should be treated with caution for some sites when 
boundary conditions are significantly different from the present day. Also, the exclusion 
of the influence of humidity and other evaporative source changes in MCIM regressions may 
be a limitation of quantifying source temperature fluctuations from deuterium excess in 
some instances. 

Citation: Lewis, S. C., A. N. LeGrande, M. Kelley, and G. A. Schmidt (2013), Modeling insights into deuterium excess as 
an indicator of water vapor source conditions, J. Geophys. Res. Atmos., 118, 243-262, doi:10.1029/2012JD017804. 


1. Introduction 

[2] The water cycle is a dynamic component of the global 
climate system and stable water isotopologs (S 18 0 and 5D, 
expressed as a deviation from Vienna Standard Mean Ocean 
Water in permil, %o, units) are valuable tracers of hydrological 
change. Isotopic variability in polar ice cores is commonly 
used to infer paleoclimatic information, such as past local 
temperature changes [e.g., North Greenland Ice Core 


All supporting information may be found in the online version of 
this article. 

'Research School of Earth Sciences, The Australian National 
University, Canberra, ACT, Australia. 

2 Now at School of Earth Sciences and ARC Centre of Excellence for 
Climate System Science, The University of Melbourne, Melbourne, 
Victoria 3010, Australia. 

3 NASA Goddard Institute for Space Studies and Center for Climate 
Systems Research, Columbia University, New York, New York, USA. 

Corresponding author: S. C. Lewis, School of Earth Sciences, The 
University of Melbourne, Parkville, VIC, 3010, Australia. (sophie.lewis@ 
unimelb.edu.au) 

©2012. American Geophysical Union. All Rights Reserved. 
2169-897X/13/2012JD017804 


Project members, 2004; EPICA Community Members, 
2006], Deuterium excess ( d) is a second-order isotopic param- 
eter linking § 18 0 and 5D and has also been used to infer 
past changes in the atmospheric water cycle at Antarctic 
and Greenland sites [e.g., Masson-Delmotte et al., 2004; 
Jouzel et al., 2007; Steffensen et al., 2008], Deuterium excess 
indicates the deviation of a given data set from the Meteoric 
Water Line [Dansgaard, 1964], which represents the linear 
relationship between § 18 0 and 5D, as shown in equation (1) 
[Craig, 1961]: 

d = 5D-$x5 ls O. (1) 

[3] During evaporation, kinetic nonequilibrium processes 
affect the relationship between oxygen and hydrogen isotopes 
and tag vapor with a d signature so that source conditions where 
vapor originates can be estimated [Vimeux et al., 2002]. Proxy 
studies adopting coisotopic approaches (d together with S 18 0 
variability) identify changes in both local temperature (Tsite) 
and in the oceanic source region temperature of precipitation 
(Tsource) [e.g., Stenni et al., 2001; Masson-Delmotte et al., 
2005a, 2005b; Stenni et al., 2010]. 
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[4] Proxy reconstructions that interpret d as an index of a 
single climatic variable (i.e., Tsource changes) adopt several 
fundamental assumptions that may not be valid under differ- 
ent boundary conditions, including that the seasonality of 
precipitation and relative humidity are largely invariant 
[Hoffmann et al., 2001; Vimeux et al., 2001]. Important 
complexities influence d and confound isotope-climate 
relationships, with field and experimental observations 
indicating that d values depend on parameters controlling 
the kinetics of phase changes during both evaporation 
and condensation, including source temperature, relative 
humidity, and wind speed [e.g., Merlivat and Jouzel, 1979; 
Jouzel and Merlivat, 1984; Cappa et al., 2003; Uemura 
et al., 2008; Luz et al., 2009). Isotope-climate relationships 
depend on the stability of the moisture source to a site, the 
air-mass distillation history, and mixing, and hence 
these may vary on different timescales and under different 
boundary conditions [Masson-Delmotte et al., 2005b; Noone, 
2008; Masson-Delmotte et al., 2011], 

[5] Various types of isotope models have helped refine 
deuterium excess interpretations. In particular, the Mixed 
Cloud Isotopic Model (MCIM) is a distillation model used 
to constrain the relationship between d and temperature 
variability [Ciais and Jouzel, 1994] and to provide a quanti- 
tative assessment of the imprint of site and source conditions 
on isotopic variability (equations (2) and (3) after Stenni 
et al. [2001]; Masson-Delmotte et al. [2005a]). Combining 
MCIM results with back-trajectory analyses of air masses 
allows atmospheric circulation and moisture source changes 
to be inferred [Helsen et al., 2006; Schlosser et al., 2008]. 
Furthermore, the incorporation of water isotope tracers into 
dynamical complex general circulation models (GCMs) 
allows atmospheric vapor transport processes to be explicitly 
resolved and climatic influences on isotopic variability at 
Greenland and Antarctic proxy sites to be examined 
[Schmidt et al., 2007; Masson-Delmotte et al., 2008; Stenni 
et al., 2010]. 

[6] Isotope-enabled GCMs fitted with tagged water vapor 
tracers are useful for analyzing moisture histories under 
varying boundary conditions [e.g., Armengaud et al., 1998; 
Delaygue et al., 2000; Werner et al., 2001]. However, 
previous tagged water-type model studies show variable source 
region changes to high-latitude sites during simulated paleotime 
slices [Werner et al., 2001; Langen and Vinther, 2009] 
and Lagrangian source approaches [Sodemann et al., 2008; 


Sodemann and Stohl, 2009) are limited in their ability to reveal 
evaporative source region variability in the past. Furthermore, 
previous GCM-based studies of the Last Glacial Maximum 
(LGM) water sources are largely based on proxy-derived sea- 
surface temperature (SST) estimates rather than simulated 
values [e.g., Delaygue et al., 2000; Werner et al., 2001]. As 
such, further research is necessary to understand past changes 
in high-latitude atmospheric circulation, moisture source region 
variability, and isotope-climate relationships. 

[7] In this study, we examine the source of changes in 
deuterium excess over select ice cores sites in Greenland 
(GRIP and NGRIP cores) and Antarctica (EDC and EDML 
cores) over a broad range of simulated climate changes to 
improve our understanding of the relationship between 
deuterium excess and climate. We utilize a novel suite of 
generalized vapor source distribution (VSD) tracers as a 
complementary set of model diagnostic tools [Lewis et al., 
2010] to examine this relationship. Region-specific numerical 
regressions (equations (2) and (3)) combined with VSD- 
enabled GCM simulations provide the opportunity to evaluate 
the robustness of deuterium excess as a proxy for source 
region temperature changes under different climatic condi- 
tions. We compare modeled d and Tsource values for ice core 
sites to synthetic d and Tsource estimates, applying the same 
numerical relationships [Stenni etal., 2001; Masson-Delmotte 
et al., 2005a] that link deuterium excess to surface conditions 
and have been used to interpret ice core records. We simulate 
mean climate conditions during preindustrial, mid-Holocene 
(6 kyr), LGM, and “water-hosing” time slices, together with 
changes occurring during simplified average El Nino- and 
La Nina-like conditions, as a broad representation of a variety 
of climatic changes. Using both simulated water isotope and 
VSD fields, we assess the skill of d as a proxy for temperature 
changes in the moisture source under various boundary 
conditions and examine the implications of seasonality for 
d interpretations. 

2. Scope of Study 

2.1. Site Descriptions 

[8] Deuterium excess is investigated at a limited number 
of ice cores sites in Greenland (GRIP and NGRIP cores) 
and Antarctica (EDC and EDML cores) (Figure 1), which 
were selected based on the availability of MCIM-based 
quantitative climatic interpretations for comparison with 



Figure 1. Map indicating location of selected Greenland (left, GRIP and NGRIP) and Antarctic 
(right, EDC and EDML) ice core sites. 
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VSD tracer results. The Greenland sites are located -320 km 
apart and differ in elevation (2917 m at NGRIP and 3200 m 
at GRIP), accumulation rate (19 cm of water equivalent per 
year at NGRIP and 23 at GRIP), and deuterium excess 
(~10.5%o at NGRIP and ~9.5%o at GRIP) [ Masson-Delmotte 
et al., 2005b], NGRIP (75°06’N, 42°19’W) is located 
downstream of GRIP (72°21’N, 38°18’W) on a common 
air-mass distillation trajectory [Masson-Delmotte et al., 
2005b]. GRIP and NGRIP provide an interesting basis for 
comparison, as dissimilarity between their isotopic ratios is 
attributed to a different mix of moisture sources to each site 
and hence regional-scale variations in advection toward 
Greenland [ Masson-Delmotte et al., 2005b]. In Antarctica, 
the EPICA Dome C (EDC; 75°06’S, 123°23’E) and 
Dronning Maud Land (EDML; 75°00’S, 0°04’E) sites also 
have different elevations (3233 m at EDC and 2780 m at 
EDML), accumulation rates (2.5 cm water equivalent per 
year at EDC and 6.4 at EDML; EPICA Community Members 
[2006]), and present-day (1960-1990 average) values of 
deuterium excess in snow (~9.3%o at EDC and ~4.7%o at 
EDML) [Stenni et al., 2010]. The air-mass trajectories to 
these two sites are influenced by different oceanic basins 
(EDC faces the Indian Ocean and EDML the Atlantic 
Ocean; Stenni et al. [2010]), and hence separate climatic 
signatures in their respective d records may be deconvolved. 

2.2. MCIM 

[9] In studies utilizing the MCIM [ Ciais and Jouzel, 
1994], the model is tuned to simulate the deuterium-to- 
temperature slope along a pathway constrained by back- 
trajectory analysis and then run with varying site and source 
climatic conditions, including surface ocean 5 18 0 (§ ls O SM ,), 
relative humidity, wind speed and site temperature, and air 
pressure [Masson-Delmotte eta!., 2004], Multiple linear regres- 
sions are then performed for site and source temperatures and 
isotopic variability over a specified region [e.g., Stenni et al., 
2001; Masson-Delmotte et al., 2004, 2005a; Stenni et al., 
2010], The inversion of simple isotopic models can then be 
used to infer site and source temperature fluctuations. 

[ 10 ] Masson-Delmotte et al. [2005a] estimate the rela- 
tive imprint of site and source conditions on isotopic 
variability over Summit, Greenland, as shown in equations 
(2-1) and (2-2): 

Ad = -0.16A T site + 0.74Ar jD „ rcB - 1.9A(5 18 0w, (2-1) 

A<5 ls O = 0.88Ar jI(c - 0.58A T source + 0.9A8 18 O.sw, (2-2) 

where A indicates the deviation from modem (core top) 
conditions. The equivalent regressions for temperature and 
isotopic variability for EDC are adapted from Stenni et al. 
[2001] as shown in equations (3-1) and (3-2): 

Ad = -0.50A T site + 1.30A T source - 2.60A<5 18 O.sw (3-1) 

A<S 18 0 = 1.0A T site - 0.61A T source + 0.95A<5 18 Osw. (3-2) 

[ 11 ] These regressions are extended to EDML, where 
coefficients of regression are reported as similar [Stenni 
et al., 2010], although it is noted that this may introduce 
minor errors into source and site estimates for this location. 
The last term (S 18 0. 5H ,) is the ocean isotopic composition. 


[ 12 ] It is not possible to calculate the full range of 
unknown source variables (temperature, relative humidity, 
and wind speed) using § 18 0 and d, so previous studies have 
applied varying coefficients describing the sensitivity of d to 
ATsite to account for uncertainties [Uemura et al., 2012]. As 
uncertainties cannot be calculated for all parameters, error 
ranges are not presented for regression coefficients, although 
these are detailed in some studies using simulations under 
different site and source temperatures [Stenni et al., 2010], 
and the order of magnitude of coefficient uncertainties is 
typically 10-20% [Masson-Delmotte et al., 2004]. Uncer- 
tainties in these coefficients include analytical uncertainty 
and uncertainties in the correction for § 18 0. sw , tuning 
parameters and the restriction of evaporative conditions to 
ATsource [Stenni et al., 2010], For the Antarctic equation 
set (3-1 and 3-1), the threshold for solid precipitation and 
the supersaturation function are tuned to correct the simula- 
tion of 1995-1996 average observed isotopic values [Stenni 
et al., 2001]. Similarly, the Greenland relationships are 
tuned to fit GRIP mean modem isotopic values, with 
adjusted parameters controlling the supersaturation of vapor 
over ice and the proportion of condensate left in the cloud 
[Masson-Delmotte et al., 2005a], These relationships also 
assume that the d value of the sea-surface water is zero 
and that a constant relationship exists between condensation 
and surface temperature. 

[ 13 ] Recent analysis demonstrates that different ranges of 
temperatures and selection of model isotopic output have 
substantial implications for the sensitivity of d to ATsite 
[Uemura et al., 2012]. Uncertainties in this sensitivity are 
attributed to changes in ATsite, ATsource, and the gradient 
between these temperatures. Given the range of <i-ATsite 
sensitivities that have been identified, the robustness of the 
isotope-temperature relationship derived from the MCIM is 
tested for a broad range of climatic changes using the isotope 
and VSD tracer modules of the GISS ModelE-R. Similarly, 
as the seasonality of high-latitude snowfall is an important 
control on isotopic variability in ice cores that cannot be 
assessed directly with ice core records [Masson-Delmotte 
et al., 2005a, 2011], these relationships are also investigated 
in the winter and summer season to examine their robustness 
to seasonal variability. 

3. GISS ModelE and Experiment Suite 
3.1. Model Description 

[14] Simulations were made using the coupled atmosphere- 
ocean GISS Mode IE -R. The version used here is the same as 
submitted to the Coupled Model Intercomparison Project 
Phase 3, with the addition of water tracers. The horizontal 
resolution is 4° latitude by 5° longitude, with 20 vertical levels 
up to 0.1 hPa in the atmosphere [Schmidt et al., 2006] and a 
13 -layer Russell ocean model of the same horizontal resolu- 
tion [Hansen et al., 2007]. Atmospheric advection uses a 
quadratic upstream scheme, with nine moments advected in 
addition to mean quantities. The ocean component is 
non-Boussinesq, mass conserving, and uses “natural” boundary 
conditions at the free surface. 

[15] Water isotope tracers ('H^O, “normal” water; 
2 H'h' 6 0 or HDO, reported as 5D;~ and 'll^O, § 18 0) are 
incorporated into the atmosphere, land surface, sea ice, and 
ocean. Water isotopes are tracked through all stages of the 
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hydrologic cycle and are adverted like water throughout the 
model, although at each phase change a fractionation is 
applied, explicitly determining equilibrium fractionation 
and with parameterizations accounting for kinetic fractiona- 
tions [Schmidt et al., 2005]. Schmidt et al. [2005] detail 
model isotopic fractionation processes, including the moti- 
vation for applying a parameterized supersaturation function 
(S = 1-0.004*T, where T is the temperature for the kinetic 
effect that occurs when condensing vapor to ice crystals in 
ice), to produce the correct relationship between 8D and 
deuterium excess in Antarctic snow. Initial oceanic isotopic 
compositions are also prescribed (8D = 8 ls O = 0%o). 

3.2. Water VSD Tracers 

[16] We apply a suite of water VSD tracers in the model, 
and atmospheric transport and condensation processes alter 
these analogously to a nonfractionating water isotope tracer 
[Kelley, 2003]. The VSD is the integrated mass of water 
vapor in each model cell and is expressed as an area integral 
of evaporative input unique to that cell. The VSD can be 
represented as a weighted sum of basis functions that are 
orthogonal to one another over the earth’s surface. This 
study uses spherical harmonics as VSD basis functions, 
which are not anchored to any particular geographic bound- 
ary and require no prior definition of regions. The surface 
source of a given member of this new suite of tracers is equal 
to the evaporation field multiplied by its associated basis 
function. The sources of water vapor are traced back through 
any cloud processes to the site of surface evaporation. The 
precipitation source distribution is a subset of the VSD, 
defined where vapor condenses to liquid. We include 144 
tracers that are resolved to wavenumber 1 1 , providing an 
effective horizontal resolution of approximately 8° x 10°. 
It should be noted that factors such as land-sea contrasts 
cause real-world precipitation source distributions to vary 
unevenly, and hence the smooth spatial VSD patterns cannot 
be interpreted precisely on small spatial scales. 

[ 17 ] In addition, we define 18 different geographic 
regions, and water evaporating from these regions is tagged 


and followed throughout the simulations. These regional 
“painted water”-type tracers (hereafter regional tracers) are 
generalized from the VSDs and used primarily to investigate 
and quantify basin-scale vapor source changes between 
simulations. The reason that painted water tracers are not 
used globally at each gridbox is that this is computationally 
inefficient and 828 distinct tracers would be required to 
obtain 8° x 10° resolution. 

3.3. Experimental Design 

[is] Additional model tracers are computationally expen- 
sive and the VSD experiments take eight times longer to 
perform than comparable experiments with only water 
isotope tracers. Thus, experiments were conducted in two 
stages. First, coupled atmosphere-ocean simulations were 
conducted for a suite of experiment to provide a broad 
range of climatic changes with which to examine isotope- 
temperature relationships (Table 1). The first of these was 
performed and run to equilibrium (100s to 1000 years) for 
the preindustrial (0k), Holocene (6k), LGM (21k), and 
water-hosing (hereafter freshwater-forced) experiments. The 
simulated surface conditions from these experiments (SST 
and sea ice) as well as other changes (greenhouse gases, 
orbital parameters, and vegetation as necessary) were used as 
boundary conditions to drive the atmosphere-only version of 
the GCM with the VSD tracers incorporated in the second 
stage of the study. 

[ 19 ] We also performed a series of simulations (El Nino- 
Southem Oscillation (ENSO+) and (ENSO-) where strong 
average idealized La Nina and El Nino surface temperature 
anomalies are imposed on the Pacific Ocean, using observa- 
tional temperatures [Reynolds and Smith, 1994] and NIN03.4 
indices over the period of 1981-2003 as surface boundary 
conditions. Monthly surface temperature anomalies associated 
with average La Nina (ENSO-) and El Nino (ENSO+) condi- 
tions were synthesized by correlating ENS03.4 indices over 
this period with observed surface air temperatures (SATs) 
and were regressed to establish mean strong anomalies that 
are applied to SAT fields as surface boundary conditions for 


Table 1 . Summary of Model Boundary Conditions for Each Time Slice 3 



C02 

CH4 

N20 


8 18 0™ 




Time (kyr BP) 

(fraction) 

(fraction) 

(fraction) 

Perihelion 

(%o) 

Ice 

Other 

VSDs 

0 

1 

1 

1 

2.85 


No 


Yes 

1 

0.98 

0.87 

0.94 

349.6 


No 


No 

2 

0.98 

0.81 

1 

332.3 


No 


No 

3 

0.97 

0.8 

0.95 

314.9 


No 


No 

4 

0.96 

0.77 

0.94 

297.4 


No 


No 

5 

0.95 

0.71 

0.95 

279.9 


No 


No 

6 

0.95 

0.71 

0.82 

263.42 


No 


Yes 

9 

0.93 

0.83 

0.89 

210.4 

0.33 

LIC 

ICE-5.2GLIC 


No 

21 (LGM) 

0.66 

0.49 

0.73 

15.51 

1 

LIC 


Yes 

21 (LGM ICE-5. 2G) 

0.66 

0.49 

0.73 

15.51 

1 

ICE-5.2G 

1 Sv freshwater (T 0°C;S 0 psu;5 ls O -30%o) 

Yes 

Hosing 

1 

1 

1 

2.85 


No 

Yes 

x 100 model years over 50-70°N N. Atlantic 

ENSO+ 

1 

1 

1 

2.85 


No 

El Nino-like average surface temperature 
anomalies 

Yes 

ENSO- 

1 

1 

1 

2.85 


No 

La Nina-like average surface temperature 
anomalies 

Yes 


a Details are provided for 0. 1, 2, 3, 4, 5, 6, 9, 21, and 21k LGM-ICE-5.2G, freshwater-forced (hosing), ENSO+, and ENSO-, including greenhouse gases 
(Indennuhle et al., 1999), perihelion in Julian days (Berger and Loutre, 1991), mean global 5 ls O.™, and ice sheet boundary conditions. 
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the Ok simulations (Supplementary Figure 1). This method 
produces surface temperature conditions that are broadly 
representative of positive and negative phases of ENSO, 
although the regression was multiplied by a factor of 3 to 
ensure that the fingerprint of this type of climatic change is 
discernible throughout the high-latitude hydrological cycle 
(Supplementary Figure 2). 

[20] The control simulation (Ok) was run with all boundary 
conditions and atmospheric composition appropriate to 
the preindustrial period (ca. 1880). Experiments were also 
performed for mean climate conditions during the mid- 
Holocene (6 kyr BP, “6k”; LeGrande and Schmidt [2009]), 
the LGM (21 kyr BP, “21k”; A. LeGrande, personal commu- 
nication), and a freshwater-forced event, as a highly idealized 
analog to a Heinrich event [Stouffer et al., 2006; Lewis et a!., 
2010], In the 6k and 21k time slices, greenhouse gas concen- 
trations were adjusted according to ice core reconstructions 
[Indermuhle et al., 1999; Brook et al., 2000; Sowers et al., 
2003] and seasonal insolation parameters [Berger and Loutre, 
1991], Surface ocean 5 ls O values are prescribed globally as 
l%o in the 21k simulation. 

[21] The ice cap used for the 21k simulation is based 
on the ICE-5. 2G reconstruction prescribed for PMIP2 simu- 
lations [ Peltier , 2004], although we employ a modified 
Laurentide Ice Sheet model (ICE-5. 2G-LIC; Licciardi 
et al. [1998]) that is considered most realistic for capturing 
LGM atmospheric circulation in the high northern latitudes 
(Supplementary Figure 1; D. Ullman, personal communica- 
tion). A second 21k simulation was run with a standard 
ICE-5. 2G [Peltier, 2004] boundary conditions as a compar- 
ison for the modified experiment. Hereafter, this experiment 
is referred to as LGM-ICE-5.2G so as to distinguish it from 
the primary simulation of the LGM. 

[22] Water-hosing simulations are forced by a freshwater 
perturbation of 1 Sv (1 Sverdrup = 10 6 m 3 /s) added over 
100 model years over 50°-70°N of the North Atlantic basin 
[Lewis et al., 2010], Each experiment was run for longer 
than 500 years in order to reach quasi-equilibrium, and mean 
climatologies over the last 100 years are present here. The 
LGM simulation was run for >800 years, although it did 
not fully reach equilibrium. Comparison of high-latitude 
surface conditions in this simulation with a complete LGM 
run confirmed the anomalies presented here are valid, 
although a secondary VSD-enabled run has not yet been 
conducted due to computational constraints. Results from 
the highly idealized freshwater-forced and ENSO-like simu- 
lations, and the LGM experiment, are not presented for 
climatic interpretations in their own right, but rather to assist 
in understanding (/-climate relationships. 

[23] For the second stage of these experiments, the atmo- 
sphere-only, VSD-enabled model was run for 6 years, with 
the final 5 years’ results reported here. A preindustrial coupled 
simulation with VSD tracers included was also conducted in 
order to test the validity of using results from atmosphere-only 
simulations, and this shows only small differences in precipi- 
tation source distributions. In addition to this suite of isotope- 
and VSD-enabled simulations, we consider an expanded 
series of Holocene time slice simulations made for 1, 2, 3, 4, 
5, and 9 kyr BP [LeGrande and Schmidt, 2009], without 
corresponding VSD simulations. These experiments allow a 
comparison of simulated isotopic fields to be made with recon- 
structed proxy changes from ice core sites. 


3.4. Definitions 

[24] During each simulation, ATsource values for each 
site for annual and seasonal (June-July-August (JJA) and 
December-January-February (DJF)) temporal averages 
were calculated as average surface temperatures within the 
simulated source region (including both land and ocean 
gridboxes, unless stated otherwise) using model VSD tracers, 
and modeled d values during these time averages are simply 
values occurring at gridboxes encompassing each ice core site. 
Corresponding “synthetic” ATsource and Ad values were 
determined from equations (2) and (3) using modeled Tsite 
and 5 18 O svv fields, as a basis for comparison with modeled 
values. Climatological conditions within VSDs are calculated 
as averages weighted by the fractional amount of vapor 
contributed to the total precipitation. 

[25] Climatic changes are reported as anomalies (A) 
relative to preindustrial (0k) simulation. Values presented 
for each ice core site (such as Ad) are those occurring within 
single gridboxes, although are coherent with changes over a 
broader area unless otherwise stated. Climatological results 
from El Nino- and La Nina-like experiments are described 
as anomalies of ENSO+ phases relative to ENSO- (AENSO). 
Southward and westward changes in mean moisture source 
locations to the sites are described by negative values, while 
northward and eastward shifts are described by positive 
values. All climatological anomalies reported are greater 
than 95% significant, given the control decadal variability 
about the 1 00-year mean. 

[26] We define recycling as water with a continental, 
rather than an oceanic, evaporative source. The mean vapor 
transport distance (TD) is calculated as the distance (km) 
between the mean location of the precipitation source distribu- 
tion and the ice core site. These values provide a minimum 
estimate of transport distances and do not consider the curved 
trajectories of air masses 

4. Evaluation of Model Performance and Modeled 
Changes 

[27] Extensive model evaluation of the GISS ModelE and 
its isotope module has been undertaken previously and 
showed that patterns of simulated preindustrial S 18 Op are 
highly correlated with observed fields [Schmidt et al., 
2005, 2006]. Large-scale simulated climate conditions and 
isotopic variability during the freshwater-forced [Lewis 
et al., 2010] and mid-Holocene [LeGrande and Schmidt, 
2009] experiments have also been assessed previously. 
Given prior model evaluation and the focus of this studying 
being the assessment of (/-temperature relationships across a 
range of different climate states, we do not comprehensively 
compare simulated climatologies against observations for 
all experiments. 

[28] The large-scale changes in simulated surface condi- 
tions across the suite of simulations are shown in Figure 2, 
and a summary of simulated site changes for the suite of 
experiments is provided in Table 2. Generally, discrepancies 
exist between simulated and observed climatological data 
in the high latitudes, with a warm bias occurring over 
Antarctica and a dry bias over Greenland resulting from 
deficiencies in the simulation of Northern Hemisphere 
(NH) storm tracks [Schmidt et al., 2006]. Modeled precipita- 
tion patterns are generally closely related to observed 
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ANN Ok SAT (°C) 13.79 



ANN A6k SAT (°C) -0.35 
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Figure 2. Annual average SAT (°C, left), relative humidity at 850 mb level (%, middle), and deuterium excess (d) in 
precipitation (%o, right) for the control simulation (0k, upper) and anomalies relative to preindustrial for 6k, 21k, and 
freshwater-forced (hosing) experiments and ENSO+ minus ENSO- (AENSO) anomalies from the coupled simulations. Solid 
gray shading indicates areas of high model topography. 


hydrological data, although simulated preindustrial snowfall 
is subdued relative to observations. These various climatic 
factors can also influence the relationship between tempera- 
ture and isotopic compositions [Masson-Delmotte et al., 
2008, 2011]. Elowever, the pattern of deuterium excess in 
precipitation is reasonably simulated, with the greatest 
discrepancies in hydrological fields occurring over the 
tropics, and isotopic values over Antarctica are generally 
enriched relative to observed [Schmidt et al., 2007]. Overall, 


we consider further analysis of simulated climates within 
the high-latitude regions to be valid, but focus on anomalies 
between experiments to account for these deficiencies. 

4.1. Comparison With Ice Core Sites 

[29] A comparison of preindustrial simulated annual average 
surface characteristics (Tsite, d, 5 1 & O p , accumulation, and 
precipitation) with observations at each of the four ice core sites 
is provided in Table 3. The simulated accumulation rate is 
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Table 2. Summary of Simulated Climatic Changes at Ice Core Sites GRIP, NGRIP, EDC, and EDML for 6k, LGM, Freshwater-Forced 
(Flosing), and ENSO-like Experiments 0 





0k 



ANNUAL 


ANNUAL 

JJA 

DJF 

A6k 

ALGM 

AHosing 

AENSO 

GRIP 

Tsite 

-29.9 

-16.5 

-39.8 

-0.1 

-12.2 

-6.2 

0.3 


PR 

0.1 

0.2 

0.1 

0.0 

-0.1 

-0.1 

0.0 


d 

11.9 

13.4 

10.1 

0.1 

-4.3 

3.9 

-0.6 


8 l8 O p 

-34.7 

-32.6 

-37.2 

0.2 

-1.1 

-4.0 

0.4 

NGRIP 

Tsite 

-30.4 

-16.5 

-39.8 

0.2 

-9.7 

-6.7 

0.5 


PR 

0.4 

0.2 

0.1 

0.0 

-0.1 

-0.1 

0.0 


d 

12.8 

13.4 

10.1 

0.2 

-4.9 

3.5 

-0.2 


5' 8 O p 

-34.9 

-32.6 

-37.2 

0.1 

-0.9 

-3.7 

0.4 

EDC 

Tsite 

-44.0 

-55.0 

-27.6 

-0.3 

-9.2 

-1.4 

0.5 


PR 

0.1 

0.1 

0.2 

0.0 

-0.1 

0.0 

0.0 


d 

6.7 

7.8 

6.0 

-0.3 

-5.3 

0.6 

0.4 


8 18 Op 

-38.2 

-40.1 

-37.1 

-0.1 

1.1 

-1.5 

0.0 

EDML 

Tsite 

-34.3 

-41.3 

-24.2 

-0.4 

-7.4 

0.4 

1.9 


PR 

0.2 

0.3 

0.1 

0.0 

-0.1 

0.0 

0.1 


d 

9.5 

9.5 

8.3 

-0.6 

-3.5 

0.8 

1.9 


8 18 O p 

-39.3 

-40.5 

-36.2 

0.3 

-0.1 

-0.3 

-0.7 


“Details are provided for SAT of site (Tsite, °C), precipitation amount (PR, rnm/day), deuterium excess in precipitation (d, %o), and 5 ls O in precipitation 
(d^Op, %o). Simulated preindustrial annual and seasonal average values are also provided. All values reported are greater than 95% significant (student’s 
t test) given the control decadal variability about the 100-year mean. 


Table 3. Comparison of Simulated Annual Average Preindustrial Site SATs (Tsite, °C), Tsource (°C), Precipitation Amount (PR, mm/ 
day), Accumulation Rate (cm Water Equivalent/Year), Deuterium Excess in Precipitation ( d , %o), and 5 18 0 in Precipitation (S 18 !},,, %o) 
With Observed Values for Each of the GRIP, NGRIP, EDC, and EDML Ice Core Sites 0 




GRIP 

NGRIP 

EDC 

EDML 

Model 

Obs 

Model 

Obs 

Model 

Obs 

Model 

Obs 

0k 

Tsite 

-29.9 

-31.6 

-30.4 

-31.7 

-44.0 

-54.5 

-34.3 

-44.6 


Tsource 

10.7 

-15-20 

9.8 

-16-21 

5.4 

-10-15 

10.7 

-5-10 


Accumulation 

4.8 

19.0 

12.3 

23.0 

2.4 

2.5 

3.2 

6.4 


PR 

0.1 

0.3 

0.4 

0.6 

0.1 

0.3 

0.2 

0.5 


d 

11.9 

10.5 

12.8 

9.5 

6.7 

9.3 

9.5 

4.7 


s 18 o„ 

-34.7 

-35.5 

-34.9 

-35.2 

-38.2 

-50.7 

-39.3 

-44.8 

21k 

Tsite 

-12.2 

-21.0 

-9.7 


-9.2 

-8.1 

-7.4 

-7.9 


Tsource 

-2.9 

-6 

-4.2 


-1.8 

-2.9 

-3.1 

-2.7 


d 

-4.3 

-6 

-4.9 


-5.3 

-2.9 

-3.5 

-2.0 


8 18 O p 

-1.1 

—5 

-0.9 


-1.1 

-5.2 

-0.1 

5.3 


“Observed values for Tsite, Tsource, d p , and 5 18 O p are provided by Stenni et al. [2010] for EDC and EDML, and isotopic values are core top averages 
(~1 .2-2.0 kyr BP). Present-day values for GRIP and NGRIP are taken from Masson-Delmotte et al. [2005b]. Simulated precipitation is compared to climate 
prediction centre merge analysis of precipitation reanalysis data [Xie and Arkin, 1996], In addition, comparisons between simulated and reconstructed 
isotope and temperature values are provided for the LGM (EDC: Stenni et al. [2004]; GRIP: Masson-Delmotte et al. [2005a]: EDML: Stenni et al. 
[2010]) and anomalies are average values over the period ~20-22 kyr BP. Modeled LGM ATsource anomalies represent average SST (rather than SAT) 
values within the simulated source region. Comparison is not given for NGRIP, as the deuterium excess record for this site has not been interpreted 
quantitatively for the LGM. 


notably lower than observed for the NH sites due to compara- 
tively smaller simulated precipitation rates. Annual average 
simulated site temperatures are similar to observed for GRIP 
and NGRIP, although the warm bias over the southern high 
latitudes is apparent for EDC and EDML. The observed 
seasonal SAT cycle at each site is captured well within the 
model (Supplementary Figure 3), with this warm bias evident 
particularly during the winter months at the Antarctic sites. 
Agreement between simulated and observed isotopic values is 
closer for GRIP and NGRIP than the Antarctic sites, where 
simulated d is underestimated at EDC and overestimated at 
EDML (Table 3). Overall, the observed intersite gradients 
(between GRIP/NGRIP and EDC/EDML) are not captured in 
the preindustrial simulation. It should also be noted that the 


Dome C region is described as localized and variable in terms 
of temperature, isotopic composition of precipitation, and 
accumulation rates [Masson-Delmotte et al., 2008], and model 
resolution is important for capturing observed isotopic gradients 
in Antarctica [Werner et al., 2001], Hence, comparison of mean 
gridbox values with observational point data may not produce a 
close agreement. Also, snowfall deposition occurring in the 
form of clear sky precipitation, which is not well resolved 
GCMs [Masson-Delmotte et al., 2008], may impact model- 
observational comparisons at EDC. 

[30] Additional comparisons with reconstructed LGM 
values show that simulated Tsite changes for EDC and 
EDML are close to the estimates from isotopic reconstruc- 
tions, although simulated Tsite and § 18 Op values at GRIP 
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Figure 3. Comparison of simulated d anomalies throughout 
modeled Holocene time slices (1, 2, 3, 4, 5, 6, and 9k) with 
reconstructed values over GRIP (a; Masson-Dehnotte et al. 
[2005b]), NGRIP (b; Masson-Delmotte et al. [2005b]), EDC, 
and EDML (c and d; Stenni et al. [2010]). Model data points 
indicate the 1 a decadal variability for the preindustrial control 
simulation as a measure of simulated natural variability. 
Estimates of analytical uncertainty of d vary, but may be up 
to ~0.64%o (Masson-Delmotte et al., 2005a). 


are subdued relative to proxy reconstructions (Table 3). 
Multiproxy approaches, including borehole reconstructions, 
suggest Tsite anomalies for the LGM in Greenland of greater 
than 20°C [Dahl-Jensen et al., 1998], while only modest 
changes of - — 10-12°C are simulated here. The model is 
able to resolve the sign of the glacial-interglacial d variations 
as observed in Greenland and Antarctic ice cores, which 
has not been captured in previous studies using simula- 
tions forced with CLIMAP SST reconstructions [ Werner 
et al., 2001]. 

[ 31 ] We also compare ModelE simulated d anomalies with 
reconstructed changes from ice core records throughout the 
Holocene (Figure 3). Over the Greenland sites, modeled 
annual average d anomalies (relative to the control simulation) 
largely resemble those reconstructed from the Holocene 
section of the GRIP record, with an increase in d from between 
9 and 6 kyr BP. The correspondence between simulated and 
reconstructed d values is lower over the Antarctic EDC site, 
with greater negative anomalies simulated in the early 
Holocene, compared with measured values [Stenni et al., 
2010], although the reconstructed d trend is simulated well 
for EDML. 

4.2. Preindustrial VSDs 

[32] For GRIP and NGRIP sites, the dominant annual 
average control source of precipitation is located in the 
western North Atlantic (Figure 4), which comprises ~72% 
of total precipitation. Regional tracers are in agreement with 
moisture sources diagnosed by reanalysis-based and concep- 
tual models [. Johnsen et al., 1989; Barlow et al., 1997; 
Sodemann et al., 2008], The precipitation source to NGRIP 
is more disperse than GRIP, and regional tracers indicate 
an increased fraction of moisture transported from distal 
sources over continental North America and the northern 
Pacific basin (by ~8%) and a decreased contribution from 
the North Atlantic (by ~14%). The increased mean vapor 
transport distance of simulated moisture to NGRIP of ~380 
km, compared to GRIP, supports previous studies indicating 
more remote precipitation sources for this site [ Masson- 
Delmotte et al., 2005b]. A greater proportion of precipitation 
is derived from the North Atlantic region southeast of 
Greenland during the boreal winter months (DIF), while 
during the summer (JJA), there is an increase in precipitation 
to Greenland sites with a Pacific Ocean evaporative source 
and from the high northern latitudes around Greenland. 
The summertime Pacific source includes a significant contri- 
bution of continentally derived recycled moisture from 
North America (~50%). 

[33] Over Antarctica, the mean source of simulated annual 
control precipitation to EDML is located at ~42°S, with the 
majority of moisture sourced from the South Atlantic and 
Southern oceans (Figure 4) There is also a contribution 
of long-distance moisture transport to this site from the 
mid-latitude Pacific in agreement with Lagrangian back 
trajectories analyses at ice core sites [Sodemann and Stohl, 
2009]. Supporting seasonal observations [Masson-Delmotte 
et al., 2008], there is a greater proportion of modeled 
moisture to EDML with a provenance in the South Atlantic 
during the austral winter months (JJA) and a greater amount 
of continentally recycled moisture during the summer 
months (DJF), when regional evaporative rates are higher. 
Over EDC, the simulated annual control precipitation 
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Figure 4. Simulated control precipitation source distributions (0k, left) and seasonal precipitation source 
region anomalies (JJA/DJF, right) for Greenland sites GRIP and NGRIP and Antarctic sites EDC and 
EDML. Solid black rectangular boxes indicate the ice core sites (end member precipitation), and dashed 
boxes indicate the mean location of the precipitation source region to each site. Precipitation source 
distributions are unitless probability density functions normalized by the maximum probability density. 


(mean location ~58°S) consists of locally sourced moisture summer-dominated snowfall regime at EDC is modeled with 
dominated by the Southern Ocean within the Indian Ocean a significant recycled component. During the austral winter 
basin and recycled moisture from continental Antarctica. A months, the main source of precipitation is comparatively 
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distal (ATD -2800 km), located within the Southern 
and southern Indian oceans (-60% from regional tracers). 
Previous studies suggest the main source of EDC precipita- 
tion is located in the mid-latitude Indian Ocean, with a more 
southward EDML moisture source region within the South 
Atlantic Ocean [Stenni et al., 2010]. However, simulated 
VSDs do not support a more southward source for EDML 
relative to EDC. Observational studies characterize EDC 
climates as highly localized [Bromwich et al., 2004], and 
distinct moisture source regimes are inferred for high- 
(EDC) and low-elevation (EDML) sites, with low-elevation 
sites more influenced by higher-latitude moisture while 
high-elevation sites are more influenced by lower-latitude 
moisture [Masson-Delmotte et al., 2011]. 

[34] Sublimation sources are an important component of 
Antarctic snowfall, with mass balance estimates suggest- 
ing sublimation effects of up to 40% in coastal locations 
[. Bromwich et al., 2004]. At EDC, however, the overesti- 
mation of recycled precipitation sources is likely, with 
back- trajectory analyses also indicating a lower-latitude 
moisture source for this site [Sodemann and Stohl, 
2009], Previous studies note that GCMs can overestimate 
the importance of recycling in the hydrological cycle, 
compared with advective moisture sources [Ruiz-Barradas 
and Nigam, 2006], which may have implications for the isoto- 
pic composition of snowfall [Masson-Delmotte et al., 2011]. 
The overestimation of recycling at EDC may result from the 
wann bias at this continental site, although it is not considered 
prohibitive to comparing synthetically derived source tempera- 
ture estimates with simulated VSD temperature values, as 
only climatological anomalies relative to the preindustrial are 
considered here. 

4.3. VSD Changes 

[35] Both VSD and regional water tracers show changes in 
source regions to the four ice core sites throughout the suite 
of simulations (Figures 5 and 6). During the mid-Holocene 
experiment, there is only a small simulated northward shift 
(relative to preindustrial) in the mean precipitation source 
to GRIP of -1° in latitude that corresponds to a -100 Ion 
decrease in vapor TD. For NGRIP, there is a similar northward 
source shift, together with a decrease in mean vapor TD, of 
-145 km during the mid-Holocene. Additionally, there is an 
increase in Pacific-sourced moisture to this site during this 
time slice (of -6%) and a concomitant decrease in the contri- 
bution of moisture originating in the Atlantic. Over Antarctica, 
a mid-Holocene southward source shift is modeled at EDML 
(—2° Alat, ATD - — 180 km), although the latitudinal change 
at EDC is comparatively smaller (—0.5° Alat, ATD —320 
km). In both high-latitude regions, there are minimal 
corresponding changes in precipitation amount during the 
mid-Holocene. 

[36] During the 21k experiment, VSD tracers indicate 
minimal mean source location changes for GRIP, with 
complex spatial VSD changes occurring (Figure 5). A larger 
mean northward source change (of -7° in latitude) is simu- 
lated at NGRIP, compared with the preindustrial simulation. 
During the LGM, both sites’ regional tracers show an 
increased fraction of North Atlantic precipitation (GRIP 
-11%; NGRIP -18%), with simulated decreases in mean 
vapor TD to these Greenland sites of up to 800 km; simu- 
lated moisture sources are more localized during the LGM. 


These changes are also associated with a decrease at both 
sites in the amount of precipitation that is greatest in 
the summertime (GRIP JJA —30%; NGRIP JJA —20%). 
Simulated SST anomalies for the LGM, relative to the 
preindustrial, are similar to synthesized proxy reconstruc- 
tions [MARGO Project Members, 2009] around the 
oceanic source regions to the northern ice core sites 
(Figure 7), although interpretations of proxy reconstruc- 
tions indicate a southward moisture source shift to GRIP 
of -5° in latitude during the LGM. Prior model studies 
assign variable source region changes to Greenland during 
the LGM under prescribed SSTs [Werner et al., 2001; 
Langen and Vinther, 2009]. The simulated reduction in 
the North American source with the addition of the 
Laurentide at the LGM (by -5%) is also observed in 
sensitivity experiments [Langen and Vinther, 2009; MARGO 
Project Members, 2009], 

[37] Over Antarctica, there are also decreases in pre- 
cipitation amount during the LGM at both sites (EDC 
ANN —60%; EDML ANN —55%). The shift in mean 
source region location for the 21k simulation is substantial 
at both sites (EDC -3° Alat, -13° Along; EDML —4° 
Alat, -2° Along) and predominantly southwestward. This 
is associated with decreases in vapor transport distances 
(EDC ATD —640 km; EDML ATD —470 km), a greater 
contribution to both sites from the Southern Ocean region, 
and a reduced mid- and low-latitude moisture component. 
Previous GCM results indicate that LGM source region 
shifts for Atlantic-influenced Antarctic regions are minimal, 
while Pacific-dominated source changes are greater [Werner 
et al., 2001], Although the strong expansion of LGM sea ice 
was highly seasonal [Gersonde et al., 2005], the southward 
LGM source shift for EDC again results in an overly signifi- 
cant contribution of sublimated snowfall (-60%) to this site 
during summer for this time slice. Considering the oceanic 
component of the EDC source only, the shift at the LGM 
is significantly westward. 

[38] Regional tracers in the freshwater-forced experiment 
reveal increases of up to -15% in the contribution of the 
North Atlantic to total annual precipitation at NGRIP and 
GRIP. This is associated with an overall decrease in precip- 
itation to these sites (- — 0. 1 mm/day), including a decrease in 
moisture advected to Greenland from the Pacific. Overall, 
the mean TD of vapor to GRIP is unchanged. Following 
freshwater injection, moisture transport distance to NGRIP 
is reduced by -134 km, compared with during the preindus- 
trial. For the EDML Antarctic ice core site, freshwater 
forcing results in a source shift close to the site (-3° Alat) 
and a corresponding decrease in mean moisture TD of 
-370 km, although there is no significant change in the over- 
all amount of precipitation to the Antarctic sites following 
freshwater injection. The Southern Ocean dominates the 
freshwater-forced source to EDML, with a decrease in 
moisture transported from the southern Atlantic. At the 
EDC site, the freshwater-forced source is largely stationary 
relative to the present day, with a small westward shift in 
mean source location by -2°. 

[39] For the ENSO-type experiments, the ENSO+ phases 
are associated with increases in annual average precipitation 
at all sites. For GRIP, there is a southward source sift 
of -1.5° and an increase in TD of -170 km. During simu- 
lated positive ENSO phases, there is a decrease in both 
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Figure 5. Simulated precipitation source distribution anomalies for GRIP (left) and NGRIP (right) sites 
for the 6k, 21k, freshwater-forced (hosing), and ENSO-like experiments. Precipitation source distributions 
are unitless probability density functions normalized by the maximum probability density. Site locations 
are indicated in Figure 4. 


Pacific-sourced vapor and moisture from the Labrador Sea to 
this site. The largest GRIP source increase for the ENSO+ 
experiment occurs from the mid-latitudes of the Atlantic. 


A similar fingerprint of ENSO phases is modeled for the 
moisture source to NGRIP. For the Antarctic sites, the 
largest ENSO simulation impacts occur for the EDML, 
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Figure 6. Simulated precipitation source distribution anomalies for EDC (left) and EDML (right) sites 
for the 6k, 21k, freshwater-forced (hosing), and ENSO-like experiments. Precipitation source distributions 
are unitless probability density functions normalized by the maximum probability density. Site locations 
are indicated in Figure 4. 


which has a stronger Pacific precipitation influence than 
EDC (~25% of total annual precipitation from the Pacific 
for EDML, ~8% for EDC from regional tracers). There is a 
shift in mean precipitation ( ~ 2 ° Alat, ~4.5° Along) and 


increase in mean vapor TD by ~550 km, associated 
with positive ENSO phases at this site. At EDC, there is an 
eastward shift in mean source by ~5° and a decrease in trans- 
port distance by -150 km. 
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A. MARGO ALGM SST(°C) 
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Figure 7. (a) Proxy-derived MARGO synthesis SST anomalies (SST, °C) (MARGO Project Members, 
2009) for the LGM and (b) equivalent SST anomalies for the ICE-5. 2G-LIC experiment using the GISS 
ModelE-R. 


4.4. Source Region Characteristics 

[40] Assessment of the accuracy of simulated average 
source region temperature changes for the core sites is 
limited, as there are few periods detailed in proxy studies 
that are characterized by significant d and source tempera- 
ture change. We consider the LGM the most useful basis 
for model-proxy data comparison of source temperatures, 
as this period is significantly different climatically from the 
preindustrial and boundary conditions are comparatively 
well known. Generally, simulated patterns of LGM SST 
anomalies are reasonably close to proxy-derived recon- 
structed values (Figure 7). At the GRIP site, the modeled 
LGM source SAT change of —1°C is close to the — 6'C 
reconstructed from coisotopic approaches [Masson-Delmotte 
et al., 2005a], Over Antarctica, the simulated LGM 
source temperature change to EDC (~7°C) is greater than 
the — 3°C suggested by proxy reconstructions [Stenni 
et al., 2010], although the modeled site temperature change 
of ^9.2°C is in good agreement with reconstructed changes 
of - — 8°C. The offset between simulated and reconstructed 


LGM source temperature changes at EDC may relate to the 
discrepancy in the proposed location of source regions 
to this site and particularly to the likely simulated overestima- 
tion of land-derived snowfall. 

5. Evaluation of d as a Tracer of Source 
Temperatures 

[41] For each ice core site (GRIP, NGRIP, EDC, and 
EDML), average modeled ATsource values determined 
using VSD tracers are compared to synthetic ATsource 
values derived from equations (2) and (3) for annual and 
seasonal (JJA, DJF) average changes for each simulation 
(Figure 8). Similarly, a comparison is made between 
modeled synthetic Ad values (Figure 8). Most modeled VSD 
average source temperature data points match synthetically 
derived equivalents well, particularly given that the uncer- 
tainty of the isotopic inversion (equations (2) and (3)) method 
is ~2-3°C on ATsource and ATsite reconstructions [Masson- 
Delmotte et al., 2005a]. 
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Figure 8. Comparison of synthetic and model VSD ATsource (panels a, c, and e) and model A d 
(panels b, d, and f) estimates for 6k, 21k, freshwater-forced (hosing), ENSO+, and ENSO- experi- 
ments for GRIP (blue), NGRIP (red), EDC (green), and EDML (grey) ice core sites. Comparisons 
are made for annual, JJA, and DJF averages. Simulated Tsource values are determined using aver- 
age values within modeled VSD regions, and model d values are determined as climatological 
values within the relevant ice core gridbox. Synthetic Tsource and d values are calculated from 
empirical relationships ( Stenni et al., 2001; Masson-Delmotte et al., 2005a) using modeled Tsite 
and S 18 0 !VV fields. The solid line indicates the 1:1 relationship between the Tsource or d changes, 
where estimates from both techniques are identical. The dashed gray line indicates the line of best 
fit for each panel. 
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5.1. Greenland Sites 

[42] For the mid-Holocene and ENSO-like experiments, a 
close agreement occurs between Ad and ATsource estimates 
from modeled and synthetic values, indicating that d is 
generally a faithful tracer of source region temperatures for 
Greenland under these boundary conditions. We identify 
the largest discrepancies in Ad and ATsource comparisons 
using these two techniques for the LGM simulation, particu- 
larly during the boreal summer. Discrepancies could indicate 
multiple climatic influences on d during this time, or 
that deficiencies exist in either MCIM-derived <i-Tsource 
relationships or the GCM experiments. 

[43] The LGM was characterized by marked topographic 
changes within the cryosphere [ Clark et al., 2009]. High- 
latitude regions are sensitive to changes in boundary condi- 
tions associated with the glacial augmentation of the ice 
sheets and sea-ice extent, such that glacial atmospheric 
circulation was significantly altered [Jouzel et al., 1982]. 
Furthermore, Greenland precipitation is a heterogeneous 
and dynamic mix of moisture sources with different distilla- 
tion histories [ Charles et al., 1994] and is responsive to 
regional perturbations in circulation patterns. During the 
LGM, the blocking of low-level atmospheric flow by topog- 
raphy (principally the Laurentide Ice Sheet) impacted polar 
front and regional cyclogenesis and advection of moisture 
into central Greenland was significantly perturbed [Werner 
et al., 2000; Langen and Vinther, 2009]. Simulated glacial 
changes in the region include an increase in eastward near- 
surface airflow around Greenland, a simulated southward 
deflection of atmospheric jets, and a strengthening of the 
polar front in response to orographic forcing (Figure 9). 
The steering influence of the Laurentide Ice Sheet on 
regional circulation also affects the moisture source and 


isotopic composition of precipitation throughout the Arctic 
region. Moisture is advected on transport pathways of signif- 
icantly different lengths and temperature profiles, which 
impacts the retention of heavy isotopes in vapor and the 
composition of end-member precipitation. Regional tracers 
indicate that during the LGM there was an annual average 
decrease in Pacific-sourced moisture (~15%) from the 
preindustrial contribution of ~24% (Figure 9). This is consis- 
tent with the reduction of the Pacific source under the 
influence of the Laurentide determined from previous model 
experiments [Charles et al., 1994]. 

[44] The largest discrepancies in ATsource comparisons 
for Greenland sites occur during the boreal summer months. 
The importance of the seasonal intermittency of precipita- 
tion for isotopic variability over high-latitude sites is 
well noted, although poorly constrained for LGM climates 
[Masson-Delmotte et al., 2005a, 2008). Glacial conditions 
inhibit winter snowfall due to both circulation changes in 
response to the Laurentide Ice Sheet and altered thermal 
oceanic gradients. Simulated dry glacial winters may result 
from a shift toward more zonal circulation patterns that 
preclude onshore advection to Greenland and reduce the 
contribution of strongly isotopically depleted precipitation 
to the annual total [Jouzel et al., 2007]. In addition, there 
is a substantial increase in the continentality of central 
Greenland during the LGM [Otto-Bliesner et al., 2007; 
MARGO Project Members, 2009). This includes a simulated 
4-8° (1-2 gridbox) latitudinal LGM expansion southward of 
sea ice in the North Atlantic in the wintertime and an annual 
average simulated increase in NH horizontal snow and sea- 
ice coverage of ~10%. Simulated LGM winter precipitation 
comprises less than 15% of the annual total to GRIP and 
NGRIP, and the seasonal cycle of precipitation is accentuated 
at this time relative to the control simulation. 


ANNA21kd(%o) -0.62 ANN A21k PRECIP (mm/day) -0.24 ANN A2 lk Jet Speed (m/s) -0.16 A21k ANN precipitation source - GREENLAND 



JJA A21k d(%o) -0.38 JJA A21k PRECIP (mm/day) -0.26 JJA A21k Jet Speed (m/s) 0.00 




DJF A2Ik d(%o) -0.86 DJF A21k PRECIP (mm/day) -0.25 

DJF A21k Jet Speed (m/s) 0.35 




> mini II !I ii m i i i i i i :|| mu. i m n m n n n m i n n n m n n n n 
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Figure 9. Simulated LGM annual (ANN) and seasonal (JJA, DJF) d (%o, left) and precipitation (mm/day, 
middle left) and jet speed (m/s, middle right) anomalies and Greenland precipitation source region anomalies 
(right). VSDs are unitless probability density functions nonnalized by the maximum probability density. The 
Greenland region (70-80°N, 35 — 45° W') encompasses both GRIP and NGRIP ice core sites. 
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[45] During summer under preindustrial conditions, the 
evaporative source to the Greenland ice core site lies 
predominantly in the Pacific Ocean and high latitudes 
around Greenland, and these sources are strongly perturbed 
by changes in the regional LGM topography (Figure 9), 
including the expanded summer sea-ice extent. For example, 
there is a ~25% increase in the contribution of the North 
Atlantic source to GRIP summer precipitation during the 
LGM. Masson-Delmotte et al. [2005a] note that the applica- 
tion of isotopic inversions without accounting for precipita- 
tion seasonality biases produces unrealistic ATsource 
estimates, and hence the significant LGM changes in the 
seasonality of both precipitation regimes and source region 
locations for the Greenland sites may result in Tsource 
discrepancies between simulated and synthetic estimates. 
In this instance, the comparison attempted here may benefit 
from seasonally specific regression coefficients. 

[46] Given the sensitivity of regional climates to topography, 
the accurate parameterization of ice-sheet orography is critical 
for accurately capturing the character of NH atmospheric circu- 
lation, precipitation [Kageyama and Valdes, 2000], sea-ice 
extent, and meridional temperature gradients [Masson-Delmotte 
et al., 2005a]. Therefore, a second comparison of ATsource 
estimates using the VSD tracers and numerical relationships 
was also made for the second LGM experiment that was run 
with standard ICE-5 .2G boundary conditions (see Table 1; 
Peltier [2004]). The correspondence of LGM ICE-5.2G 
Tsource changes for the two techniques is significantly lower 
than for the modified Laurentide equivalent (average Greenland 
5.2G-LIC Tsource R 2 = 0.7; 5.2G Tsource R 2 = 0.2). This 
difference highlights the importance of topography on regional 
hydrology and indicates that ATsource and d estimate offsets 
may be an artefact of poor model ice-sheet representation. It is 
possible that in this topographically sensitive experiment, the 
model resolution and prescribed LGM boundary conditions 
utilized here may be inadequate to resolve topography critical 
for regional climatologies, and a higher-resolution model incor- 
porating ice-sheet dynamics, such as horizontal advection, 
may improve the simulation of regional LGM climatologies 
[Carlson et al., 2008] and the comparison of ATsource 
estimates. However, it should be noted that the correct sign of 
LGM present-day d anomalies is captured in the simulations 
utilized here, whereas in previous atmospheric LGM simula- 
tions forced LGM forced with CLIMAP SSTs, erroneous 
d anomalies in Greenland occurred [e.g., Werner et al., 2001]. 

[47] Overall, comparison of Tsource estimates at these 
Greenland sites suggests that deuterium excess may poorly 
record source temperature conditions in some instances by 
the proposed empirical scalings (equations (2) and (3); Stenni 
et al. [2001]; Masson-Delmotte et al. [2005a]), particularly 
where boundary conditions are substantially different. As 
certain aspects of simulated LGM hydrology over Greenland 
(e.g., A d p over GRIP and VSD-derived ATsource) do corre- 
spond closely with reconstructed proxy records, it may be that 
simplifications in numerical J-Tsource primarily contribute to 
offsets between the ATsource estimates for the LGM. For 
example, the pronounced LGM seasonality changes may 
affect the applicability of numerical regressions (equations 
(2) and (3)) in the manner attempted here, as these are valid 
for individual precipitation events, but do not take into account 
seasonality changes in either snowfall or the isotopic composi- 
tion of precipitation [Jouzel et al., 2007]. The assumption of 


temporal invariance in seasonality is known to result in 
erroneous site temperature estimates from isotopic profiles 
from Greenland and seasonal corrections have been 
required previously [Masson-Delmotte et al., 2005b]. 
Although the isotope inversion models have been applied 
at glacial-interglacial timescales, their validity for boundary 
conditions very different from the present remains to be 
demonstrated, and in these cases, unique inversions may be 
necessary [Uemura et al., 2012]. 

5.2. Antarctic Sites 

[48] For EDML, Tsource and d changes are similar for 
modeled and synthetic estimates (Figure 8), indicating that 
deuterium excess may provide a robust proxy for moisture 
source temperatures in different climatic states. Conversely, 
larger disparities occur at EDC, the higher-elevation Antarctic 
site, particularly for the LGM and summer freshwater-forced 
comparisons. Synthetic and simulated EDC ATsource esti- 
mates are similar, although corresponding Ad comparisons 
are comparatively poorer. Comparisons for the EDC site 
may be susceptible to poor simulation of snowfall regimes 
due to a simulated warm continental bias and an overestima- 
tion of moisture recycling from the continent. Simulated 
preindustrial d p is, for example, enriched relative to observed 
values in this region (Table 2). In this region, d is particularly 
topographically dependent, and model-observation mismatches 
may result from inaccuracies in representations of both cloud 
microphysics and water-vapor advection {Masson-Delmotte 
et al., 2008). Discrepancies in Tsource and d estimates suggest 
that this site is indeed susceptible to local or regional-scale 
model shortcomings. 

[49] Although Ad offsets are seasonally robust for EDC 
and largely ubiquitous across the suite of simulations, the 
largest discrepancies occur at the colder end of the climate 
spectrum (i.e., the LGM and freshwater-forcing experi- 
ments). As such, inadequacies in simulated climatological 
regimes at this site (such as the noted warm bias) may be 
exacerbated during these extreme cold times. For example, d 
variability depends on microphysical processes prevailing in 
clouds, including kinetic fractionation, during the processes 
of condensate formation from atmospheric vapor [Jouzel and 
Merlivat, 1984], which are variable under differing climate 
states. For snow formation occurring in a supersaturated 
environment, such as over ice, a kinetic fractionation effect 
transpires similar to that existing during evaporation ( Jouzel 
et al., 1 982, 1991) and depends on highly local supersaturation 
conditions when condensing water vapor to ice [Schmidt et al., 
2005]. For a given deuterium excess in vapor, the effect of ice 
condensation from atmospheric vapor is to lower d values 
in snow, although this is also temperature dependent {Luz 
et al., 2009). The impact of condensation processes during 
snow formation is significant, particularly at the cold end of 
the climatic spectrum, and on the order of several per million 
{Luz et al., 2009). 

[50] Over Antarctica, where AT source estimates correspond 
more closely with each other than over Greenland, the impact 
of cold climate conditions in Antarctica is considerably more 
subdued. There are minimal changes in the seasonality of 
simulated LGM precipitation over Antarctica, which largely 
maintains its geographical symmetry during glacial condi- 
tions. Prior model studies indicate the effect of seasonality is 
of minor regional importance on glacial precipitation regimes 
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[Jouzel et al., 2003; Stenni etal., 2010). Previously, Antarctic 
d variability has been considered amenable to simpler interpre- 
tations than for Greenland [Jouzel et al., 2003]. In this study, 
the qualitative correspondence of source temperature and 
d values is closer for EDML than at other sites and this 
isotope-temperature relationship appears robust to changes in 
boundary conditions. 

6. Discussion 

[ 51 ] In this study, we examined changes in the SAT of the 
evaporative source to various locations, while reconstruc- 
tions using ice core coisotopic variability generally interpret 
fluctuations in terms of source SST changes. Although SSTs 
are considered generally equivalent to SATs near the ocean 
surface in observational isotope analyses [Uemura et al., 
2008], significant SST-SAT offsets (up to ~4.5°C during 
the 21k and freshwater-forced experiments) occur in 
modeled source regions using VSD tracers. Studies based 
on climate reanalyses show that although there is a close 
SST-SAT correlation in most regions, this relationship 
degrades significantly in the high latitudes, where the 
difference in air-sea temperatures is influenced by multiple 
variables, such as wind speed, surface heat fluxes, and net 
solar radiation [Cavan, 1992; Kara et al., 2007]. In cases 
where source air-sea temperature differences are large, for 
example, changes in humidity and latent heat fluxes at 
the evaporative surface may also be significant, with larger 
air-sea temperature contrasts causing higher heat fluxes 
[Song and Yu, 2012], 

[52] Indeed, there is a degraded correspondence of 
J-ATsource values when SST, rather than SAT, source 
averages are used for comparison between model and 
synthetic estimates (Table 4). Average source SAT values, 
when compared with SSTs, integrate a larger combination 
of influences on the isotopic composition of vapor during 
evaporative processes. Additionally, recycled moisture 
originating from land, rather than oceanic, sources is 
an important component of precipitation to proxy sites, 
including those in the high latitudes, and also influences 
the composition of water vapor. Disregarding this mois- 
ture component by investigating source SSTs, rather 
than SATs, neglects a significant contribution to air-mass 
histories. 

[53] Source region characteristics that control evaporation 
kinetics, such as relative humidity and wind speed, also 
provide important climatic influences on isotopic composi- 
tion [Merlivat and Jouzel, 1979]. For example, d correlates 
with both relative humidity at the ocean surface and with 
source temperatures [Uemura et al., 2008] In this study, 


we model significant changes in relative humidity between 
simulations and a seasonally dependent correspondence 
between d and relative humidity changes at the Greenland 
sites (Table 4). There are average relative humidity changes 
of ~25% within the Greenland area precipitation source 
region during the LGM (Figure 9) and smaller simulated 
LGM anomalies (~3-9%) for the Antarctic sites. Notably, 
simulated Greenland source humidity differences during 
the LGM (Figure 9) are an order of magnitude greater than 
source relative humidity changes accounted for in Tsource-r/ 
equations (2) and (3) and may introduce errors in numerical- 
based reconstructions. The exclusion of the influence of 
humidity changes in MCIM regressions may be a significant 
limitation of quantifying source temperature changes from 
isotopic profiles. 

[54] Additionally, the geographical location of the evapo- 
rative source relative to the site is an important influence on 
the isotopic composition of end-member precipitation, as the 
proximity of the moisture source impacts the degree of air- 
mass mixing and rainout incurred en route to the proxy site 
[Charles et al., 1994], Evaporative source regions are found 
to be dynamic and there are significant simulated transport 
distance reductions from source to site between the prein- 
dustrial and LGM simulations (of up to ~800 km), which 
may influence the isotopic composition of end-member 
precipitation, as moisture is advected on different pathways 
and sustains a varying degree of replenishment by other 
moisture sources. Deuterium excess is dependent on advec- 
tion height — lower d values transpire in near-surface flow 
originating from localized high-latitude evaporative sources 
and higher-moisture transport pathways reflect more distal 
transport [Stenni et al., 2001; Masson-Delmotte et al., 2008]. 

[55] Generally, deuterium excess variability is associated 
with changes in multiple climate parameters and the relative 
importance of these is both regionally dependent and 
susceptible to changes in prescribed boundary conditions. 
We find that simulated Greenland hydrological regimes are 
sensitive to varying boundary conditions. Given the similar- 
ity of simulated and synthetic ATsource estimates for the 6k, 
freshwater-forcing, and ENSO-like simulations, it may be 
that beyond a certain climatic threshold (i.e., during the 
LGM), multiple climatic influences should ideally be 
considered in interpreting d variability. The simplifications 
incorporated into the numerical inversions (equations 
(2) and (3)), particularly that the seasonality of precipitation 
is temporally static and that the source relative humidity is 
largely invariant, may introduce a bias into isotope-climate 
relationships. These considerations may ultimately result in 
the offset between synthetic and modeled ATsource values 
over Greenland during the LGM. 


Table 4. Correlation Coefficients (R) for GRIP, NGRIP, EDC, and EDML Site d and Simulated Source Region Average SATs and SSTs 
and Relative Humidity for Annual (ANN) and Seasonal (JJA, DJF) Values Correlation Coefficients Were Determined From Modeled 
Values Across the Suite of Simulations 


Site 

d- 

Tsource (SAT) R 

d- 

Tsource (SST) R 

d- Relative humidity R 

ANN 

JJA 

DJF 

ANN 

JJA 

DJF 

ANN 

JJA 

DJF 

GRIP 

0.4 

0.8 

0.0 

0.2 

0.3 

0.0 

0.0 

0.1 

-0.2 

NGRIP 

0.6 

0.8 

0.2 

0.2 

0.3 

0.1 

-0.1 

0.1 

-0.3 

EDC 

0.6 

0.5 

0.7 

0.3 

0.4 

0.2 

0.1 

0.1 

0.1 

EDML 

0.7 

0.6 

0.7 

0.6 

0.6 

0.5 

0.3 

0.3 

0.1 
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6.1. Limitations of the Study 

[ 56 ] There are limitations to this study that preclude the 
deuterium excess from being fully assessed as a moisture 
source tracer under different climatic states. First, the simpli- 
fications inherent in the numerical equations (2) and (3) 
mean that several important climatic parameters are not 
explicitly represented. The regression coefficients are 
regionally specific, but are determined over a range of 
climatic states and are not boundary condition specific 
[Uemura et al., 2012]. Changes in boundary conditions 
impact the permanence of the evaporative source, precipita- 
tion seasonality, and elevation corrections, together with 
evaporation and condensation processes. Flence, refined 
equations incorporating further climatic variables (such as 
relative humidity), or providing boundary condition or season- 
ally specific coefficients, may reduce the discrepancies 
between source temperature estimates determined syntheti- 
cally from MCIM regressions and those simulated here using 
average values within the VSD. This would be particularly 
enlightening for Greenland model-proxy comparisons for 
changes occurring in the LGM experiment. 

[57] Second, although GCMs are useful tools for examin- 
ing isotope-climate relationships because they capture 
dynamical complexity beyond Rayleigh-type processes, 
they are, by necessity, only as good as the accuracy of 
parameterizations (such as kinetic fractionation processes) 
relative to the real world. The parameterizations of the 
MCIM and GISS models are different and so conducting 
MCIM simulations using parameterizations may help 
resolve the discrepancy in source temperature estimates 
identified here. Furthermore, model biases, such as poorly 
resolved inversion layers, excessive accumulation, and 
warm continental biases, may result in erroneous simula- 
tions of polar climates [Schmidt et al., 2007]. Orography in 
GCMs is strongly dependent on both the model resolution 
and discretization scheme ( Masson-Delmotte et al., 2006). 
Hence, the model resolution utilized here may be inadequate 
to resolve critical topography for regional climatologies and 
a higher-resolution model may improve the simulation of 
precipitation regimes at sites such as EDC, which is highly 
sensitive to local and regional-scale topography. An increase 
in horizontal and vertical model resolution would likely 
yield improved d performance over Antarctica [Werner 
et al., 2011], In addition, the combined use of 17 0-excess 
with deuterium excess may help disentangle source effects 
on isotopic signals [Winkler et al., 2012], 

7. Conclusions 

[58] Using the water isotope-enabled GISS ModelE-R 
and a novel suite of vapor source distribution tracers, we 
examine the relationship between deuterium excess variability 
at high-latitude ice core sites and the climatic conditions 
prevailing within their moisture source regions. Average 
modeled source temperature values for various Greenland and 
Antarctic sites are compared to synthetic values calculated from 
numerical isotope-temperature relationships (equations (2) and 
(3); Stenni et al. [2001]; Masson-Delmotte et al. [2005a]) for 
a suite of simulations encompassing a range of mean climatic 
states. These include highly idealized experiments that are 
included here to assess the rigor of (/-temperature relationships 
under varying boundary conditions. Similar comparisons are 


made with simulated d changes over ice core site gridboxes 
and synthetic d values. 

[59] There is a general correspondence between modeled 
and numerically derived source values for most sites, particu- 
larly during the mid-Holocene and ENSO-type climate simu- 
lations, indicating that in some simulations and regions, d is 
a faithful tracer of temperature variability within a site’s evap- 
orative source using MCIM-derived inversions (equations 
(2) and (3)). Over Greenland sites GRIP and NGRIP, there 
are larger discrepancies between Tsource estimates occurring 
during the LGM experiment, particularly in the boreal 
summer. Simulated LGM climates are characterized by 
substantial changes in regional circulation under the influence 
of expanded sea ice and the Laurentide Ice Sheet. The quanti- 
fication of source temperature fluctuations may be affected by 
changes in precipitation intermittency and evaporative mois- 
ture source origins during the LGM. The disparity of modeled 
and synthetic estimates of Greenland Tsource changes during 
the LGM demonstrates the complexity of (/^-Tsource relation- 
ships, particularly under vastly different boundary conditions 
where specific regression coefficients may be necessary 
[Uemura et al., 2012], There are also large discrepancies in 
comparisons of synthetic and modeled d changes at the EDC 
site in Antarctica, which are attributed to regional scale 
deficiencies in modeled climatologies. Here significant dispa- 
rities exist between simulated preindustrial climates and local 
observations, indicating that this high-elevation, topographi- 
cally complex area requires high spatial resolution modeling 
to capture realistic climatologies [Masson-Delmotte et al., 
2011; Werner et al., 2011], Conversely, at the lower-elevation 
EDML site, deuterium excess was found to be a faithful tracer 
of evaporative source temperatures by the numerical isotope- 
temperature relationships utilized and is a generally robust 
proxy throughout a range of different climatic states. At this 
location, coisotopic approaches are invaluable for reconstruct- 
ing past changes in the high-latitude hydrological cycle under 
a broad range of conditions. 

[60] There are multiple well-recognized climatic influences 
on d variability, including those that drive moisture evapora- 
tion, condensation, and transport. We suggest that the closer 
correspondence of modeled and synthetic source SATs, 
compared with SST averages, results from the influence of 
evaporation kinetics on the isotopic composition of water 
vapor. Changes in relative humidity and wind speed are 
important climatic influences on isotopic composition, and 
significant changes in relative humidity are modeled here in 
different experiments. The exclusion of the influence of hu- 
midity changes in MCIM regressions may be a limitation of 
quantifying source temperature changes from isotopic profiles 
where large changes in humidity in the source region occur, 
such as during the LGM. Similarly, continentally derived 
moisture provides a significant contribution to overall precip- 
itation and its composition, even for high-latitude proxy sites. 
The apparent importance of a range of climatic controls on 
evaporation supports results from previous model, experimen- 
tal, and observational-based studies that highlight the contri- 
bution of relative humidity at the ocean surface and the 
strength of winds to deuterium excess fluctuations. 

[61] The climatic changes simulated here in a series of 
idealized experiments using a relatively coarse horizontal 
resolution GCM, provide a useful complement to observa- 
tional and isotopic model-based studies for understanding 
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deuterium excess as a climate proxy and high-latitude paleo- 
hydrological processes more generally. Deuterium excess is 
an integrated tracer of past hydrological changes and records 
a suite of climatic changes. The numerical equations resulting 
from simple isotopic models are a simplification of the climate 
system, and many variables that affect the isotopic composi- 
tion are not explicitly parameterized. As a proxy, d is found 
to be amenable to simplification at some sites and for 
some types of climatic change. However, variability in other 
controls of evaporation kinetics, such humidity and wind 
speed, or in precipitation seasonality and moisture source 
region, are also important isotopic influences and these 
warrant further investigation. Quantitative interpretations of 
deuterium excess may require caution for some sites at times 
characterized by significantly different boundary conditions 
than the present, particularly during the LGM. 
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